Non-adiabatic effects of superconductor silane under high pressure 
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Investigations of non-adiabatic effects by including vertex corrections in the standard Eliashberg 
theory show that high phonon frequency is unfavorable to superconductivity in regime of strong 
vertex correction. This means that it is hard to find high-transition-temperature superconductors in 
the compounds with light elements if the non-adiabatic effects are strong. The interplay interaction 
between non-adiabatic effect and Coulomb interaction makes the transition temperature of silane 
superconductor not so high as predicted by the standard Eliashberg theory. 
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The recently found superconductor silane (S1H4) at 
high-pressure [l[ partially realizes the theoretical predic- 
tions of high-temperature superconductors of hydrogen 
dominant metallic alloy @ and metallic hydrogen at high 
pressure However, the transition temperature (T c ) is 
significantly lower than the values predicted in previously 
theoretical calculations 0-0] • What is the underlaying 
mechanism of the lower T c compared with the higher 
T c predicted in strong-coupling theory. With increas- 
ing pressure the electronic structures of silane will evolve 
from molecule-like energy-levels to electron-bands when 
gaseous silane transforms into crystalline silane. The ef- 
fective band-widths of conducting electrons of metallic 
silane are narrow so that the non-adiabatic effects will 
be prominent. It's very desirable study the influences of 
non-adiabatic effects on T c . 

The non-adiabatic effects are partially included in the 
vertex corrections in the theory of electron-phonon in- 
teraction. The Eliashberg theory [8l-[Toj combined with 
the vertex corrections had been widely studied by per- 
turbation method pdl4l3j . We have performed a full 
parameter-space search based on Eliashberg theory with 
and without vertex corrections In this paper, 

we study the influences of the vertex corrections on T c 
of superconductor silane at high pressure using reliable 
Eliashberg functions a 2 F(ui) obtained from the calcula- 
tions of linear-response theory Q. Our results indicate 
that the vertex correction can efficiently suppress T c ap- 
proaching to the values found in experiments. Addition- 
ally, we find that high phonon frequency is unfavorable 
to high-T c if the strong vertex corrections are included. 
This means that it is hard to realize high-T c supercon- 
ductors in silane and other hydrogen-rich materials. 

We have generalized the equation of energ y g ap in 
Ref.[ll| by including the Coulomb interaction 13]. The 
calculations of vertex corrections are greatly simplified 
under isotropic approximation. The electron-phonon in- 
teractions are included in the vertex corrections only by 



the functions A n defined as A„ = 2 J °° dua 2 F \v)v / \v 2 + 
cj 2 ). When temperature is very close to T c , the 
energy-gap equation including the leading vertex cor- 
rection from electron-phonon interaction is written as 
'" 1 A n r / 1 ui n > I = with the kernel matrix 
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where the definitions of H n , A nn i, B nn >, C nn ' and a n 
are presented in Ref.flU Il3j . In order to calculate T c , 
the matrix K nn > is symmetrized with the same man- 
ner as m Ref.[l0|. The E B is the effective band-width 
and the Coulomb pseudo-potential fi* is defined as //* = 
Mo/(l + MoM-Es/Oo)), where ^ = N(0)U, N(0) the 
density of state at Fermi energy Ep, U the Coulomb in- 
teraction between electrons and f2o characteristic energy 
of typical phonon correlated to superconductivity. Gen- 
erally, the ratio f2o /Eb takes as the parameter to measure 
vertex correction and the larger Qq/Eb is correspond- 
ing to the stronger vertex correction. In this context, 
the vertex corrections are controlled by Eb, the smaller 
Eb (or larger Qq/Eb) for stronger vertex corrections. 
With development of method of electronic-structure cal- 
culation in solid materials, the parameters of electron- 
phonon interaction can be calculated using density func- 
tional theory combined with perturbing linear-response 
theory [l5|, [l^]. The Elaishberg functions a 2 F(uj) ob- 
tained by linear response theory are easily merged in the 
formulae used in this context. 

In the standard Eliashberg theory, T c will infinitely 
increase with phonon frequency or energy. It's widely 
expected that high-temperature superconductors should 
be found in compounds containing light elements such as 
hydrogen-rich material silane and metal hydrogen at high 
pressure. We will see that the non-adiabatic effects set 
bounds to infinitely increasing T c . In order to obtain very 
general results, at first step, we apply the simple phonon 
spectrum and Eliashberg function a 2 F(uj) used in Ref.0]. 
The Figd shows how T c changes with fip and Eb- We 
can see that for large band- width Eb~5 eV (weak vertex 
correction), T c increases monotonously with flp just as in 
the standard Eliashberg theory. However for very strong 
vertex correction with small band-width Eb=1 eV, T c is 
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in Ref [4] in the calculations of linear-response theory. 
The structures of crystal silane at high-pressures haven't 
been completely defined. The metallic P63 structure 
had been found in experiment [lj. The theoretical lay- 
ered structure with Cmca space-group symmetry is more 
stable from 60-200 GPa other stable structures at 
higher pressures had been already reported [l7j . The 
hydrogen-rich superconductor silane at higher pressure 
has very high phonon energies coming from the vibra- 
tions of hydrogen atoms. The effective phonon energies 
(u))i n for silane at high pressures distribute from 50 meV 
to 75 meV dependent on pressures. The parameters A of 
electron-phonon interaction are about 1.17, 0.62, 0.75 for 
the Cmca structure and 0.84, 0.87, 1.1 for the C2/c(2D) 
structure respectively when pressure increases from 70 
GPa to 150 GPa and to 200 GPa. The T c calculated using 
Allen-McMillan formula are distributed from 20 K to 80 
K Q that are greatly larger than experimental T c < 17K. 



FIG. 1: The T c on Qp — Eb parameter plane. The two solid 
arrows A and B show the changes of T c with flp at fixed Eb 
(A) near 1 eV and (B)5 eV respectively. The dashed arrow 
shows the direction enhanced the vertex corrections. 




FIG. 2: The dependence of T c on Eb and fi* , based on the 
Eliashberg functions a 2 F(ui) obtained from linear-response 
theory at different pressures for the Cmca structure (a,b) 
and C2/c(2D) structure (c,d) respectively 3. The (a) and 
(c) include vertex corrections (VC), and there are no vertex 
corrections (NVC) in (b) and (d). 



non- monotonously dependent on fip. When 51p is larger 
than a threshold value, T c will decrease with ftp. This 
means that, if the vertex corrections are included, high 
phonon energy is unfavorable to superconductivity. 

We preform the calculations of the vertex corrections 
of real superconducting material silane containing the 
lightest element: hydrogen. The lowest order vertex 
correction can significantly reduce T c when energy of 
phonon Hp is larger than 80 meV shown in FigfTJ We 
adopt the Eliashberg functions that had been reported 
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FIG. 3: The T c on [i*-Eb plane using the Eliashberg func- 
tions o?F(io) obtained from linear-response theory j4|. The 
structure has C2/c(2D) space group symmetry at pressure 200 
GPa. The shade areas are regions with T c <17 K. The arrow 
shows the directions of increasing the vertex correction. 

As reported in Ref.[ij], the pressure-dependence T c of 
the Cmca structure is different from the C2/c(2D) struc- 
ture. The T c of the Cmca structure decreases when pres- 
sure increases from 70 GPa to 150 GPa, and slightly 
increases with increasing pressure further to 200 GPa. 
However for the C2/c(2D) structure, T c monotonously 
increases with pressures from 70 GPa to 200 GPa. We 
first present the results of the Cmca structure in detail 
shown in FigfJ]Ja,b). If the vertex corrections are not in- 
cluded and ;U*=0.1, the values of T c are close to 20 K 
at pressures 150 GPa and 200 GPa, which are smaller 
than those obtained from McMillan formula. So our 
results indicate that for very broad distributed phonon 
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spectrum, the standard strong-coupling calculations are 
needed to accurately calculate T c beyond simple McMil- 
lan and Allen T c formulae. 

At low pressure with larger inter-molecule distances, 
the crystalline silane is an insulator with narrow energy 
bands corresponding to discrete energy levels of a single 
molecule. The effective band- widths will increase with 
increasing pressures and overlaps of molecular orbits of 
different molecules. If the vertex corrections are included 
in our calculations with the effective band- width Eb of 
conducting electrons from 0.5 eV to 5 eV, the T c at 150 
GPa and 200 GPa are significantly reduced to the value 
smaller than 17 K shown in FigfSJa). The effective band- 
widths obtained from the density functional calculations 
increase from very small values at low pressure to about 
1.0-3.0 eV at 100-200 GPa for the C2/c(2D) and Cmca 
structures [4|]. In another case, the effects of Coulomb 
interaction are more significant because the T c decreases 
to very small values with increasing /i* at 150 GPa and 
200 GPa shown in Figf^b). If we want to know what 
is more important for the depression of T c we should 
know real value of /z* at 150 GPa and 200 GPa. At 70 
GPa, the small band- widths with Eb<1 eV guarantee 
the T c lower than 17 K. We find that, if without vertex 
correction, even larger fi* (=0.4) can not suppress T c to 
value smaller than 17 K at 70 GPa. 

The effects of vertex correction of the C2/c(2D) struc- 
ture in FigfSfc) are similar to the Cmca structure in 
FiglUa), however the influences of the Coulomb interac- 
tion in Figdjd) are different from the Cmca in FigfJIh). 
The C2/c(2D) structure is less stable than the Cmca 
structure. The stronger electron-phonon interaction in- 
duces its structure more unstable at high pressures from 
150 GPa to 200 GPa. Theoretical values of T c increase 
with pressures. This is because the enhancements of T c 
induced by the increase of A are larger than the depres- 
sions of T c due to the vertex correction. For the C2/c(2D) 



structure, the Coulomb interaction depresses T c to values 
which arc still larger than 17 K shown in Figf2][b). So 
Coulomb interaction individually can not explain why T c 
smaller than 17 K in high-pressure experiments. The ef- 
fective band-widths Eb of conducting electrons in the 
C2/c(2D) structure are larger than 1 eV at 150 GPa 
and 200 GPa [H so that the vertex correction individ- 
ually can not explain low T c at 100 GPa and 200 GPa 
as well as shown in FigfJ^c). The Figj3] illustrates the T c 
for the C2/c(2D) structure on \i*-Eb plane at pressure 
200 GPa. The shading region in the figure shows the 
effective region with parameters \i* and Eb that can ex- 
plain the Tc smaller than 17 K, which is only small area 
of parameter-space with ^i*>0.1 and Eb<2.0 eV. Espe- 
cially, the hatched region shows the parameter-space to 
explain why experimental T c is smaller 3 K at 200 GPa. 
So our results indicate that the interplay of vertex cor- 
rection and Coulomb interaction can explain the low T c 
of the C2/c(2D) structure. Our calculations can not fully 
explain pressure-dependent T c in Ref [l[, especially, the 
higher T c =17 K near 100-120 GPa because of the lack of 
information of structural changes. 

In summary, we study the effects of the vertex correc- 
tion on T c for superconductor of silane at high pressure. 
Our results indicate that the non-adiabatic effects are 
the barrier to prevent T c from increasing infinitely with 
phonon frequency. This means that it is hard to real- 
ize high T c to home temperature in high pressure rich- 
hydrogen materials. Our results also show that the inter- 
play interaction between Coulomb interaction and vertex 
correction is essentially to explain the larger differences 
between theory and experiments. 
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